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We propose a dipole-force linear waveguide which con- 
fines neutral atoms up to A/2 above a microfabricated single- 
mode dielectric optical guide. The optical guide carries far 
blue-detuned light in the horizontally-polarized TE mode and 
far red-detuned light in the vertically-polarized TM mode, 
with both modes close to optical cut-off. A trapping min- 
imum in the transverse plane is formed above the optical 
guide due to the differing evanescent decay lengths of the 
two modes. This design allows manufacture of mechanically 
stable atom-optical elements on a substrate. We calculate the 
full vector bound modes for an arbitrary guide shape using 
two-dimensional non-uniform finite elements in the frequency- 
domain, allowing us to optimize atom waveguide properties. 
We find that a rectangular optical guide of 0.8 /im by 0.2 /im 
carrying 6mW of total laser power (detuning ±15 nm about 
the D2 line) gives a trap depth of 200 fiK for cesium atoms 
(piF — 0), transverse oscillation frequencies of f x = 40 kHz 
and fy — 160 kHz, collection area ~ 1 /im 2 and coherence 
time of 9 ms. We discuss the effects of non-zero tuf, the Dl 
line, surface interactions, heating rate, the substrate refrac- 
tive index, and the limits on waveguide bending radius. 
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There has been much recent progress in the trapping 
and cooling of neutral atoms, opening up new areas of 
ultra-low energy and matter- wave physics 13]. Waveg- 
uides for such atoms are of great interest for atom op- 
tics, atom interferometery, and atom lithography. Mul- 
timode atom waveguides act as incoherent atom pipes 
that could trap atoms, transport them along complicated 
paths or between different environments, or deliver highly 
localized atom beams to a surface. Single-mode waveg- 
uides (or multimode guides populated only by atoms in 
the transverse ground-state) could be used for coherent 
atom optics and interferometry HfjJ, as well as a tool 
for one-dimensional physics such as boson-fcrmion dual- 
ity |(|^2|,^3| and low-dimensional Bose-Einstein conden- 
sation effects |24| ]. 

The optical dipole-force has long been used to trap and 
manipulate atoms jy as well as dielectric particles [||. 
The available intensity of lasers has allowed a multitude 
of such atom traps in the far-detuned regime, giving very 
low decoherence and heating rates, and storage times on 
the order of seconds JsJ . 

Evanescent light waves have been popular in many 



atom mirrors, traps and guides [ 1C| , |11| ] since they can pro- 
vide potentials with high spatial gradients (decay lengths 
~ A/2-7T where A is the optical wavelength), and use rigid 
dielectric structures (prisms, fibers) to define the poten- 



tial shape. For example, there has been a series of repul- 
sive (blue-detuned) evanescent-wave (EW) traps which 
rely on gravity to provide the counteracting force p"2fl and 
recent experiments have shown that hollow optical fibers 
can guide atoms confined within the hollow core using a 
repulsive evanescent field guided by the fiber (ll]]2(]] . 

The idea of using an EW to provide both attractive 
and repulsive forces is due to Ovchinnikov et at |ll|, who 
proposed the use of two colors (i.e. red and blue detun- 
ings) and differing evanescent decay lengths to achieve a 
trap with the potential minimum a distance ~ A from a 
prism surface. Until now, this design has been restricted 
to planar traps (weak confinement in the other two di- 
mensions) . 

In this paper we discuss a two-color trap based on the 
EW fields above a single-mode, submicron optical 'chan- 
nel' waveguide. The trap provides tight confinement in 
two dimensions and allows free de Broglie wave propaga- 
tion in the third, forming an atomic waveguide that could 
transport atoms a between A/4 and A/2 above the opti- 
cal guide surface. Our proposal is to utilize the differing 
vertical evanescent decay lengths of the two polarizations 
carried in the single- mode optical guide (see Figure [l|). 
The physical origin of this decay length difference is the 
fact that the TM mode is closer to optical cut-off than 
the TE mode at the same frequency. 

Our proposal is reminiscent of some existing resonant 
enhancement schemes for EW mirrors (demonstrated 
with surface plasmons [jl4| and dielectric waveguides p5|) 
but with a radical change from a planar geometry to a lin- 
ear, the mechanism for exciting the guide, and the simul- 
taneous guiding of a second frequency of opposite detun- 
ing. It also shares the feature of two guided colors with an 
atom trap proposal using microsphere whispering-gallery 
modes [fL6| . 

Our design has many desirable experimental features: 
1) very little optical power is required to obtain large 
trapping intensities since the optical bound mode has 
very small cross sectional area (~ 0.3 nm 2 ), 2) the op- 
tical field is non- divergent, so can be maintained over 
distances orders of magnitude further than diffraction- 
limited propagation in free space allows, 3) the trapping 
potential is well-known, mechanically stable, and insen- 
sitive to experimental parameters other than the optical 
powers, since it is defined by single-mode intensity dis- 
tributions fixed relative to a substrate, 4) fabrication of 
arrays of closely spaced atom waveguides is possible [[3l| , 
for parallel lithography or measurement, creating "on- 
chip" integrated atom-optical elements, 5) the atoms are 
exposed providing additional optical and physical access 
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(a feature not shared with hollow-fiber designs), and 6) 
the velocity of the atoms along the direction of the waveg- 
uide could be controlled by standing waves in the light 
carried by the optical guide j2M. 

Compared to a Zeeman-effect magnetic trap for neutral 
atoms, far-detuned optical dipole-force traps can have 
comparable trapping times, but typically an order of 
magnitude less depth and transverse mode spacings than 
recent magnetic traps [Q^- However, in microfabricated 
applications the stray magnetic fields decay as a power 
law with distance, whereas evanescent light fields decay 
exponentially (ignoring for now any scattering into free 
space caused by optical defects). We believe this could 
give guided optical traps a distinct advantage in terms of 
achievable density of independent atom-optical elements 
on a single substrate. 

Also, optical traps have the advantage that there is no 
significant loss mechanism which can remove atoms from 
the trap (assuming the thermal energy is much less than 
the trap depth): spontaneous events cause a small heat- 
ing rate, and non-adiabatic changes in mp can change the 
optical potential but not the fact that the atom remains 
trapped. This is to be contrasted with a non-adiabatic 
spin-flip event in a magnetic trap, which results in loss of 
the atom. This makes optical waveguides particularly at- 
tractive for incoherent transport, when the loss of coher- 
ence due to the spontaneous events is unimportant. Fi- 
nally, optical manipulation has the advantage over mag- 
netic manipulation in terms of high possible switching 
speeds. 

This paper is organized as follows. In Section | we de- 
scribe the dipole potential, the exponential approxima- 
tion for the EW fields, and the mechanism for the differ- 
ence in decay length. We show how we optimized the op- 
tical guide dimensions, in the case of a rectangular guide 
on a substrate of unity refractive index (for rap = 0), 
and discuss some design objectives and implementation 
issues. In Section [n] we give simulated results for ce- 
sium atoms: trap depth, coherence time, transverse mode 
spacing and Q factor, and spontaneous heating rate. We 
also show how depth and coherence time are generally 
limited by only two parameters (the detuning and the 
normalized decay length difference). We study both the 
case of a substrate refractive index of unity, and in Sec- 
tion IIC the more realistic index of 1.32. We describe 
the numer ical electromagnetic finite element technique 
in Section III, including the accuracy achieved. Section 
[lV| is an investigation of two potential causes of loss or 
decoherence of atoms, namely interactions with the di- 
electric surface and bending of the waveguide. Finally 
m Section we conclude and give some future prospects 
for this proposal. 



I. TRAP CONCEPT 



A. Theory of the light potential 

An atom in a near-resonant light field of frequency ui 
experiences both a conservative force (due to stimulated 
photon exchange) and a dissipative force (due to spon- 
taneous photon emission) Jl7],|l8| . The conservative force 
is the gradient of a spatially-dependent potential £/di P (r) 
which can be viewed as the time-averaged induced dipole 
interaction energy (proportional to the real part of the 
classical polarizability) in the electric field, or equally 
well as the 'light shift' (that is, energy level shift due to 
the ac Stark effect) of the atomic ground state 

We assume that we apply a monochromatic light field 
of detuning A = lo — loq to an alkali atom (with the 
ns ^ np transition resonant at loq), in the far-detuned 
regime (A greater than the excited state hyperfine split- 
ting, but much less than loq) and the low saturation 
regime (i7 ra bi <C A, where the Rabi flopping rate O ra bi is 
defined || by ftil ra bi = I^Eq, the dipole matrix element 
multiplied by the electric field amplitude). The dipole 
potential has both a scalar and a magnetic part: 

TiT V 

Udip (r)=p s — -s(r) 

+ A i Bohr5nuc(£,S', J,i, F)m F \H(r)\ , (1.1) 

where T is the spontaneous decay rate, s(r) is the 
spatially-dependent saturation parameter, and the po- 
tential is taken to be much less than the ground state 
hyperfine splitting. Only the magnetic part is affected 
by mp , which is defined as the projection of the total 
angular momentum F on the direction of the local effec- 
tive magnetic field H(r). The constant g nuc is the nuclear 
Lande g-factor appropriate for the F of the ground state. 
The scalar potential is identical to the case of a two-level 
atom, apart from the strength factor (3 S which is i for 
detuning from the Dl line, | for the D2 [||. 

It is important to be precise with the definition of the 
saturation parameter. We write 



S (r) 



— -p2 



(1.2) 



where Eq sat = 2/ sat /eoc (in the MKSA system) is the 
squared electric field amplitude in a plane wave of in- 
tensity /sat Q ■ This avoids the ambiguities present with 
the usual definition s(r) = /(r)/J sa t in the case of a gen- 
eral monochromatic light field (composed of an arbitrary 
coherent sum of travelling and evanescent waves), and 
emphasizes that it is the local electric field alone that 
causes the dipole potential. 

The effective magnetic field has a strength and direc- 
tion given by the circularly polarized component of the 
electric field [||, which can be written thus: 

„, x a ^ r jft[Eg(r)] x 3[E (r)] 
^BohrH(r) = /3 m — — ^ , (1.3) 
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where the strength factor (3 m is — | for Dl, | for D2, 
and we take the physical electric field (with amplitude 
Eo(r) = |Eo(r)|) to be the real part of a complex field 
E(r,i) = Eo(r) exp(— iujt). The re ason for the 'dummy' 
constant /iBohr is aesthetic, so that (1.1) can be expressed 
in a standard magnetic form. Note that for nonzero mp 
the magnetic and scalar contributions to the potential 
are of similar order, if the fields have significant circular 
polarization (this will be true for our trapping fields, for 
the reason that the optical guide is close to cut-off). 

The fact that J7di P (r) has its sign controlled by the sign 
of the detuning allows both attractive (red-detuned) and 
repulsive (blue-detuned) potentials to be created. The 
potential scales as // A but the spontaneous emission rate 
scales as //A 2 ; from this follows the well-known result 
that, if coherence time is an important factor, it is best 
to be far off-resonance and use high intensities in order 
to achieve the desired trap depth yj^. 

For simplicity, in this paper we will restrict our further 
analysis and simulations to mp = 0, although our initial 
calculations suggest that the effect of the magnetic part 
of our potential when trapping in other mp states will 
not pose major problems (assuming the spin axis adia- 
batically follows the H(r) field direction), and can even 
be used to our advantage by increasing the depth and the 
transverse oscillation frequency uj x in the case tuf > 
pq| . Also, we will consider the effect of only a single 
resonance (choosing D2 because it has a larger (3 S than 
Dl), which is a valid approximation when the detunings 
from this resonance are much less than the alkali atom 
fine structure splitting Af s . Even when it becomes ad- 
vantageous to use a large detuning of the order of Af s , it 
is possible to cancel the effect of the other resonance by 
a simple shift in the two detunings (as we will see at the 
end of Section II A). 



If we now have two light fields of differing frequency, 
the atomic potentials add [|l^,^3), as long as we as- 
sume that the timescale of atomic motion is much slower 
than the beating period (that is, the inverse of the fre- 
quency difference). In our case, atomic motion occurs 
at ~ 10 5 Hz and our light field frequency difference is 
~ 10 15 Hz, so this assumption is valid. Choosing equal 
but opposite detunings ±A about the D2 line, the trap- 
ping potential for mp — is 

2 KT r 

C / di P (r) = - — — [sbiue(r) - s re d(r)] , (1.4) 

written in terms of the saturation parameters for the two 
colors. The spontaneous scattering rate j?J is a factor 
A/r smaller than t/ c jip( r )/^ (in fact this relation applies 
for any mp state), but it depends on the sum of the 
saturation parameters rather than the difference. It is 
also spatially dependent, and has the form 

2 r / r \ 2 

r sca tt(r) = - - ( J [sbiuc(r) + s rcd (r)] . (1.5) 



B. Design of the light fields 

Our basic task is to create intense evanescent light 
fields with a potential minimum sufficiently far from a di- 
electric surface to make the surface interaction potential 
and h eating mechanisms negligible (discussed in Section 
IV A ). The main difficulty arises because the evanescent 
fields have a typical exponential decay length ~ A/27T, so 
if we are to have a trap of useful depth, we are restricted 
to keep it within roughly A of the surface (less than a 
micron) . 

A potential minimum in one dimension can be obtained 
using a blue (repulsive) light field of higher intensity at 
the dielectric surface than the red (attractive) light field, 
and ensuring the decay lengths obey L re d > ^biuo, giving 
a potential of the form 



U dip (y) = A bluc e-y/ L ^ - A lcA e-y/ L ^ 



(1.6) 



This gives a repulsive force at short range, which be- 
comes attractive at long range (see Figure |l]b), and is 
the scheme for the planar trap of Ovchinnikov et al. 
]l3[ . A large amount of insight into our proposed trap 
can be gleaned from this simple one-dimensional model 
(which we call the exponential approximation), because 
the squared electric fields above the guide will turn out to 
approximate exponential forms in the vertical direction 
quite closely. 

If we define a normalized decay length difference = 
(i rc d — £biue)/£biue, then we can give two reasons why 
increasing «l is a vital design objective. Firstly it is 
easy to show that for small ap. the deepest available trap 
depth (found by optimizing the ratio of surface intensi- 
ties A re d/Abiue) scales as qj,. Secondly, a larger is 
beneficial for trap coherence, (giving a smaller sponta- 
neous decay rate at a given trap depth and detuning), 
because t he su m of the intensities can be kept lower (see 
equa tion (1.5)) for a given intensity difference (equation 



II A 



£L4)). We will quantify this latter connection in Section 



Our two key differences from the proposal of Ovchin- 
nikov et al. are as follows. Firstly, we create a non-zero 
cth by using two orthogonally-polarized bound modes of 
a dielectric slab guide, which have different evanescent 
decay lengths at the same frequency 26 1 . This contrasts 
with Ovchinnikov et al. who suggest varying the decay 
lengths by varying the reflection angles from the inside 
surface of a glass prism. Secondly horizontal confine- 
ment is achieved by limiting the width of the slab guide 
to approximately A (forming what is called a channel 
guide [p7|), which automatically creates a maximum in 
each light intensity field in the horizontal direction. This 
results in a tight horizontal confinement in the atomic 
potential, of similar size to the vertical confinement, and 
is something very hard to achieve in a prism geometry. 

A schematic of our design is shown in Figure [l]a. The 
optical guide height H and width W are kept small 
enough to guarantee that there are exactly two bound 
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modes, differing in polarization but not in nodal struc- 
ture (in optical terminology this is called single-mode): 
Ejn has an electric field predominantly in the x direction 
P5fl , and is to be excited by blue-detuned laser light, and 
E\ x has electric field predominantly in the y direction 
and is to be excited by red-detuned laser light. We can 
see why their vertical decay lengths differ by consider- 
ing the case of the slab (i.e. taking the width W — > oo), 
where these modes are simply the slab TE and TM modes 
respectively. For both these slab modes the purely trans- 
verse field obeys the differential equation 



= fc o 2 [(n«) 2 -n(y) 2 ]0 ; 



(1.7) 



where 4> = E x (H x ) for the TE (TM) mode, the eigen- 
value rigg = kg / k is the effective refractive index for the 
i th mode (k z being the wavenumber in the propagation 
direction and k the free space wavenumber), and n(y) 
is the spatially-dependent refractive index |27j. (This is 
equivalent to a one-dimensional quantum problem in the 
direction normal to the slab, in a potential —k 2 n(y) 2 with 
h 2 /2m — 1). However the boundary conditions on the 
slab surfaces differ for the two mode types: <fi is always 
continuous, but d<fi/dy is continuous for TE as opposed 
to n~ 2 d<fi/dy continuous for TM. This asymmetry exists 
because the permittivity e — n 2 varies in space but the 
permeability \i is assumed to be constant. This disconti- 
nuity in the gradient for the TM mode forces it to have 
a lower n Q g than that of the TE mode, which means it 
is less tightly bound so has a longer evanescent decay 
length. This effect becomes more pronounced as the slab 
index increases or as optical cut-off is approached (which 
happens when n c ff decreases until it reaches n s and the 
mode becomes unbound). This tendency is preserved 
even as the width is decreased to only a few times the 
height, as in our scheme. 

No analytic solution exists for the general rectangular 
guide, so we used the finite element method discussed in 
Section 111 to solve for the bound mode n e ff values and 
fields as a function of guide dimensions. Figure shows 
the resulting 'cut-off curves', that is, contours of constant 
n e ff in the parameter space (W,H). In this example we 
chose a guide index n g — 1.56 (typical for a polymer 
dielectric) and, as a preliminary substrate index 

n s = 1. 

The single-mode region, in which we wish to remain, 
is bounded below by the E\ x and Ef t curves and above 
by the Sfi curve. Note that, as in any dielectric guiding 
structure uniform in the z axis, the lowest two modes 
(E\ 1 and Efi in our case) never truly reach cut-off, 
rather, they approach it exponentially as the guide cross- 
section is shrunk to zero. For this reason, we chose 
the practical definition of cut-off for these modes to be 
n G ff = 1.05, which corresponds to only about 20% of the 
power being carried inside the guide. In contrast, higher 
modes do have true cut-offs [^8 27 (this distinction is il- 
lustrated by the dispersion curves of Figure ||), and for 



the E\ x mode our (numerically limited) contour choice of 
n c ff — 1.02 falls very close to the true cut-off curve. 

Using the numerically calculated electric 
field strengths of the Ef t and Ef x modes, we found the 
red and blue guided laser powers which gave the deep- 
est trap, subject to the constraint of fixed total power 
(keeping the detuning constant [^9|). We also imposed 
the restriction that the zero of trapping potential come 
no closer than 100 nm along the vertical line x = 0, which 
kept the trap minimum a reasonable distance from the 
surface (see Section IV A). Performing this optimization 
over a region of the parameter space covering the single- 
mode region gave a contour plot of maximum achievable 
depth for a given total power, shown within the rectangle 
overlayed on Figure |^. This depth increases from negli- 
gible values in the top left to the largest depths in the 
lower right, indicating that choosing W and H to be in 
this latter corner of the single-mode region is best for 
depth. The depth shows very little variation with W in 
this corner, rather it is clear that varying H to stay within 
our definition of the single-mode region has become the 
limiting factor on achievable depth. We indicate a prac- 
tical choice of W = 0.97 A and H = 0.25 A as a small 
marker on Figure ^|. Example trapping potentials shapes 
possible with these parameters are shown in Figure |3|; we 
discuss their properties in Section 

In Figure |l|a the direct excitation of the optical guide 
by the two laser beams is shown only schematically. In a 
realistic experimental setup this coupling into the guide 
would happen on the order of a centimeter from the atom 
guiding region, and could involve tapered or Bragg cou- 
plers |2^] from beams or from other fibers. At this dis- 
tance we estimate that isotropic stray light due to an 
insertion loss of 0.5 would have 8 orders of magnitude 
less intensity than the EW fields in the guiding region. 
Assuming the light is coherent, this limits the fractional 
modulation of the guiding potential to 10 -4 . More im- 
provements are possible, including the use of absorbing 
shields, bending the guide through large angles away 
from the original coupling direction, and reducing the 
coherence length. 



C. Discussion of optical cut-off and substrate choice 

We have shown in Figure ^| that parameters optimized 
for trap depth are near optical cut-off. It is worth gather- 
ing together the physical reasons for this. Firstly, on gen- 
eral geometric grounds, the typical available intensities in 
a guide scale inversely with the effective cross-sectional 
area of the bound mode, which, far from cut-off (that 
is, when n g — n e ff <C n g — n s ) follows very closely the 
cross-sectional area of the guide. This makes it favorable 
to shrink guiding structures to areas less than a square 
wavelength, where they generally become single-mode. 
Secondly, once we're in the single-mode regime, as we 
approach cut-off the mode power is carried increasingly 
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outside the guide, increasing the ratio of surface intensity 
to guide center intensity. Thirdly, the evanescent decay 
length in the vacuum is longer as we approach cut-off (re- 
call that in the case of the slab, this is exactly expressed 
by L^ 1 = 2fco\/ 71 cff — 1 @l> w ^ ere the factor of 2 arises 
because we are considering intensity decay length rather 
than amplitude). In the special case of n s = 1, the decay 
length diverges to infinity as we approach cut-off. Finally, 
the ratio «l becomes larger as we approach cut-off, with 
corresponding beneficial effects on depth and coherence 
(due to in creas ing the 'goodness factor' we will introduce 
in Section II A). 

Unfortunately, these purely theoretical reasons for ap- 
proaching cut-off are in opposition to more practical ones. 
The closer to cut-off a guide is, the more sensitive it is 
to manufacturing variations in cross-section: in our case 
this will be predominantly a sensitivity to H . The re- 
sult is that small variations in H cause large variations 
in mode size, or, at worst, complete cut-off. If the mode 
size change is rapid (nonadiabatic) along the z axis, (for 
instance if this change is due to surface roughness or re- 
fractive index inhomogeneities) the resulting mismatches 
will be a source of scattering of the guided power. Any 
coherent scattering back down the guide will set up pe- 
riodic modulations of the light field over long distances. 
(One way to reduce the distance over which coherent ad- 
dition is possible is to use very broad line-width light 
sources, which would dramatically reduce modulations 
due to both guided and stray scattered light). The conse- 
quence for the atoms would be a z-dependent trap depth 
and shape, and this could lead to partial reflection or even 
localization of the matter waves. In general we expect 
scattering to limit how close to cut-off we can reliably 
operate. 

With regard to the substrate, further practical issues 
arise. In the above cut-off calculation we chose the 
simplest case of n s = 1, corresponding to a guide sur- 
rounded by vacuum. A real substrate with n s > 1 has 
the unfortunate effect of limiting the propagation con- 
stant k z of strictly bound modes to be larger than the 
freely propagating wavevector in the substrate; in other 
words, 7i c ff > n s must hold or the light field will rapidly 
tunnel into the 'attractive potential' of the substrate. 
This in turn limits the decay lengths and oil that can 
be achieved. The trap properties quoted in the abstract 
and in Sections [I A and |IIB| paper rely on very low n e g 



values (1.07 for the E\ x mode, 1.18 for Ef-^) for the rea- 
son that a low n e g is the only way to create long evanes- 
cent decay lengths in the vacuum. (This is equivalent to 
Ovchinnikov et al. choosing reflecti on a ngles very close 
to critical [jl3).) However, in Section |II C| we present pre- 
liminary results for a substrate of sodium fluoride (the 
lowest-index common optical mineral, at n = 1.32), and 
do not believe the substrate alters the basic feasibility of 
our waveguide. 

For completeness, here we list some other possible ap- 
proaches to the substrate issue. 1) Use an aerogel sub- 
strate, which can have exceptionally low refractive indices 



and low loss (films of several fim thickness with indices 
of about 1.1 can be produced Q). 2) Use a dielectric 
multilayer substrate with an effective index of unity or 
less (very low loss multilayer mirrors |}0| with effective 
indices less than unity can be created). 3) Investigate 
if there exist guide shapes which have sufficiently small 
tunnelling rate into a conventional substrate that the fact 
that the modes are not strictly bound becomes irrele- 
vant (for instance, a wedge shape with the smallest face 
in contact with the substrate). 4) Unsupported guiding 
structures could be produced over short distances [ j3l"| . 
Finally, it is important to note that the idea of replacing 
the substrate by a metallic reflective layer is not practical 
because they are too lossy. 

Ultimately, the best values of W and H, the best guide 
cross-sectional shape, and the substrate choice will de- 
pend on many of the above factors and is an area for 
further research. 



II. TRAP PROPERTIES 

In the bulk of this Section we will examine the atomic 
waveguide properties for light nearly resonant with the 
D2 line of cesium, using an optical guide of index 1.56 
of the dimensions W ~ 0.97 A and H = 0.25 A from Sec- 
tion IB, and a substrate of unity index. The saturation 
intensity for cesium is 11.2 W/m 2 jf), 



and its resonant 
wavelength of 852 nm requires that the physical guide 
size is 0.83 /im by 0.21 fim. (At the end of the Section 
we present preliminary calculations for a n s — 1.32 sub- 
strate and a different guide, and discuss how the atom 
waveguide properties are changed). 

Given the guide, we are free to choose three experi- 
mental parameters, namely the optical powers carried in 
the two modes, and the detuning A (assumed to be sym- 
metric, that is, to be of equal magnitude for red and blue 
beams, because little advantage can be gained with an 
unsymmetric detuning). The first two of these can use- 
fully be reexpressed as total power P to t = -Prod + P>iue, 
and the power ratio p = Pcd/Pbiuc- The trap shape 
will be affected by p alone: we show the trapping poten- 
tial shapes achievable at the two practical extremes of 
p = 0.4 and p = 0.2 in Figure ||, where we have chosen 
Ptot and A to give identical trap depths and coherence 
times. Smaller p values cause the trap minimum to move 
further from the surface (a distinct advantage), to be 
less "bean" shaped (i.e. to have smaller cubic deviations 
from a 2D harmonic oscillator), and to cause a slight in- 
crease in collection area. It is possible to achieve a trap 
minimum as distant as yo = 0.52 A from the surface when 
p = 0.2. The only disadvantage to implementing these 
smaller p values is that a higher P tot is required to achieve 
the same trap depth and coherence time (for instance a 
factor of 7.5 increase is required as we take p from 0.4 to 
0.2). This can be quantified within the exponential ap- 
proximation, and it can be found that the total power re- 
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quired to maintain a given depth and coherence time with 
a fixed trap geometry scales as P to t ~ (1 + p)/p 1+1 ^ aL ■ 

If we were purely interested in maximizing trap depth 
at a given P to t and detuning, it would be best to make 
p as large as possible, however if we take p much larger 
than 0.4 the trap is brought so close that the corners of 
the "bean" shape touch the dielectric surface (see Figure 
||, upper plot) and we will lose effective collection area 
due to sticking of atoms onto this surface. 



A. Depth, coherence time, and Q factor 

We may ask what trade-offs are necessary between trap 
depth and coherence time. It turns out that, within the 
exponential approximation (|Oj), this is elegantly quan- 
tifyable. We can define a 'goodness factor' 
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Sbiuc(ro) - Sred(r ) _ r 
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where the trap minimum position ro is at (x = 0, y = 
and the second equality is verified by substitution of (1.4) 
and (L5), and defining U max = |£/dip(i"o)| and t co1i 



r~ att (r ). We use this latter definition because we are 
interested in the coherence time of atoms spending time 
close to the trap minimum (which will c erta inly be true 
for the transverse ground state.) Using (1.6) to solve for 
yo and evaluate the 'goodness factor', it turns out that 
the factor is independent of either laser power (i.e. of 
giving 



G 



Lied ~ -C'blu 
Lied + Lblv 
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Combining (2.1) and (|2.2[) gives 



^max^coh 
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(2.3) 



fixing the product of achievable depth and coherence time 
as a constant multiple of the detuning. This is a remark- 
able result since it shows that increasing ai is really the 
only objective in the field design of two-color EW traps. 
We can write this in units more convenient for cesium 
trap design, thus, 



t/max . Z^L = (611.2) " L 
/iK ms 2 + qj, 

-<1B)-W, 
nm 



nm 



(2.4) 



where the value ah — 0.47 ± 0.02 (taken from best-fit 
exponentials to the numerically-found squared electric 
fields for the guide dimensions of Section [I]| with n s = 1) 
has been substituted to give the the final form. This de- 
sign expression does not give the P iot required to reach 
a desired balance between C/ max and r co h, however, the 
total laser power is usually in the mW range, several 



orders of magnitude less than in most free-space trap de- 
signs. For instance, with P t ot = 20 mW, p — 0.4 and 
A = ±15 nm we could generate a trap of 2mK depth 
with the relatively short coherence time of 0.9 ms. The 
transverse oscillation frequencies in this trap would be 
u> x /2ir = 116 kHz and uj y /2n — 490 kHz (the field shapes 
fix this ratio at about 1:4), giving an atomic mode spacing 
due to the x motion of 5.6 ^K, roughly twice the cesium 
MOT temperature, and a characteristic ground-state size 
of 26 nm by 12 nm. 

For coherent guiding, we can define a more physically 
meaningful figure of merit, Q = cjj_T co h, which tells us 
the typical number of coherent transverse oscillations we 
can expect multiplied by 2tt (i.e. it is the Q-factor of 
the transverse oscillations). We should choose uj± = u> x 
since this is the smaller of the transverse frequencies in 
our case. For Q » 1 the transverse atomic modes will 
be well resolved, and o ur g uide can be a useful interfer- 
ometric device. Using (2^3), in conjunction with the fact 
that when the trapping potential shape is fixed then lj x 
is proportional to the square root of the depth, tells us 
that for a given trap and detuning, Q oc l/co x . For a 
higher Q we should choose smaller transverse oscillation 
frequencies, that is, shallower traps. For example, the 
2 mK trap discussed above has Q w 650, but if we re- 
duce it to a 20 fiK trap of the same A (by changing p or 
the laser powers), the Q is 10 times larger. Increasing A 
would allow even higher Q to be rea lize d. 

The dependence on detuning in (2.4) is another way 
of expressing the advantages already known about us- 
ing far off-resonant beams [0,0] ■ However, our single- 
resonance approximation will break down if the detuning 
is too large: we have (somewhat arbitrarily) chosen a 
detuning limit of 15 nm, as compared to Af s = 43 nm 
for cesium. At this limit, the additional dipole potential 
created due to the detunings from the Dl line is very sig- 
nificant. However, by removing the detuning symmetry 
(changing A + from +15 nm to +12.07 nm and A_ from 
-15 nm to -17.14 nm), the desired D2 single-resonance ap- 
proximation potential is recovered in the true physical 
situation of both resonances present. (These required 
shifts, which are of order A 2 /Af s , can easily be found us- 
ing the expression for the sum of dipole potentials from 
the two lines.) An additional necessity for our limit is the 
fact that any larger detunings start to demand separate 
bound-mode calculations for the two colors, a treatment 
we reserve for the future. This detuning limit in turn lim- 
its the depths, coherence times and Q-factors we quote 
here, but we anticipate similar future EW atom waveg- 
uide designs which explore the region |A| > Af s (or even 
|A| ~ luq), and achieve much better coherence. 



B. Other properties 

We estimate the collection area of the trap as the cross- 
sectional region within which the potential is deeper than 
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the typical cesium MOT energy ^bTmot- For our exam- 
ple 100 /iK traps of Figure || this area is about 1/mi 2 . 
However, it is not possible to do much better than this 
with our design: if one tries to increase the area by 
increasing the overall trap depth, the ^bTmot contour 
touches the substrate, indicating that atoms at this en- 
ergy can reach the substrate surface, where they will 
stick, limiting the effective collection area. 

To investigate the lifetime of atoms transported in- 
coherently (the multi-mode regime), we can calculate 
the heating rate along similar lines as Grimm and Wei- 
demiiller g. We start with their equation (23) which 
gives the rate of change of the average of the total en- 
ergy of atomic motion E — Ek in + E pot 



as 



E = k B T R T s 



(2.5) 



Tr being the recoil temperature, and use the assumption 
that in an equilibrated 3D trap E^ in = ^ksT. Since there 
is harmonic motion in two directions but free motion in 
the third, the virial theorem gives us E po t = l-Ekin- Com- 
bining this with r sca tt = U max T /GhA from (|2.l| ) gives 
the heating rate 



T 



5G A 



-T, 



(2.6) 



which is of the order of one recoil temperature per coher- 
ence time. For our 100 /iK depth trap at A = 15 nm the 
rate is 4.4 /iKs -1 , implying that storage and transport 
for many seconds is possible. For simplicity, we have ig- 
nored the fact that there may be distinct longitudinal and 
transverse temperatures which do not equilibrate over the 
trapping timescales. 



C. Effect of a realistic substrate 

In this section we present calculations, performed using 
the method of Section III, for a practical substrate choice 
of sodium fluoride (the lowest refractive index common 
mineral, with n s = 1.32 at a wavelength of 852 nm), and 
investigate how this changes the atom waveguide prop- 
erties from those presented above. We increased n g to 



1.7 (dense flint glass, e.g. BaSF type) in order to provide 
sufficient index difference from the substrate. 

Fixing the width at W = 1.00 A, we found that a height 
H = 0.34 A gave the largest oil of 0.20±0.01, and allowed 
both modes to be sufficiently far from cut-off (greater 
than half the power being carried inside the guide for 
both modes). The result is a goodness factor G which 
is approximately half that of the n s = 1 case, with a 
corresponding halving of the achievable product of depth 
and coherence time according to (2.4), and doubling of 
t he heating rate at a given U max and A according to 
fl2.6p . The shorter decay lengths of 56 nm and 68 nm 
(compared to 93 nm and 137 nm for n s = 1) cause the 
typical trapping distance yo to be reduced by a factor of 
roughly 1.8. 



We found that in order to reproduce the depth of 
100 /iK and yo = 0.24 /mi of the upper trap of Figure 
H (with A unchanged) we needed P to t = 22 mW, giving 
"?coh = 9 ms. The large power increase over the lrhW 
required for n s = 1 is explained by the fact that this 
yo is now towards the upper limit practically achievable 
rather than the lower. (If yo is instead scaled in pro- 
portion to the new decay lengths, the required increase 
in P to t is only a factor 1.7). In this example, we find 
the transverse oscillation frequencies have increased to 
u> x /2ir = 81kHz and uj y /2n = 202 kHz, compared to the 
original u} x /2tt = 26 kHz and Lo y /2n — 109 kHz. The in- 
crease in ujy is explained entirely by the shorter decay 
lengths, and the increase in uj x (by a factor of over 3) is 
attributed to tighter optical mode shapes. It is clear that 
this latter effect outweighs the decrease in r co h , implying 
that the inclusion of the substrate has actually increased 
Q by 50%. 

In summary, the effects of including a realistic sub- 
strate limit the maximum trapping distance yo that can 
be achieved (because of the reduction in decay lengths) , 
lower the goodness factor, increase the heating rate and 
the required optical power, but also increase the oscilla- 
tion frequencies. For our substrate choice, each of these 
changes was approximately a factor of 2, and we believe 
that they do not alter the basic practicality of imple- 
menting our proposed waveguide. 



III. NUMERICAL SOLUTION OF THE LIGHT 
FIELDS 

The detailed electric field distribution is very impor- 
tant in calculating the trapping potential above the 
waveguide. An approximation to the form of the fields in 
the y direction is given by the analytically-known solu- 
tion for the slab waveguide, but to get more accuracy and 
knowledge of the full potential shape in the x-y plane, we 
used a full-vector finite element calculation. 

The technique represents the electric and magnetic 
fields as simple piecewise functions over many "elements" 
(regions) subdividing a slice through the guide and sur- 
rounding media in the xy plane, therefore by a finite num- 
ber of degrees of freedom. Each element has a dielectric 
constant associated with it, allowing arbitrary stepwise 
refractive index distributions in the xy plane to be mod- 
elled. Maxwell's equations for propagating solutions of 
the form exp(i/c z z — iu>t) are reduced to a generalized 
sparse eigenvalue equation with fc 2 as the eigenvalue and 
the bound mode field distributions as the eigenvectors 
|p5| . Specifically, we used the technique of Fernandez 
and Lu |53| , with H x and H y as the field degrees of free- 
dom, for simplicity using first-order (bilinear) functions 
to represent these fields over a non-uniform but separable 
rectangular grid of elements. This required a generaliza- 
tion of the Fernandez and Lu implementation, and care- 
ful consideration of their line-integral terms (which are 
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non-standard for a finite element formulation) |36|]. Of 
the many available finite element approaches to dielectric 
waveguide mode solving, this frequency-domain method 
was chosen for its absence of 'spurious modes', its ability 
to handle index step discontinuities, its small number of 
required degrees of freedom and its matrix sparsity J34|. 

Rather than emulating a radiative boundary condition 
(a notoriously hard task usually requiring an iterative 
procedure due to the k z dependence), we enclosed the 
problem in a large, perfectly-conducting box of sufficient 
size that the bound mode evanescent fields were negli- 
gible on its walls, making the nature of the boundary 
condition irrelevant. However, the average level spac- 
ing of the unbound modes (the 'continuum') decreases 
with increasing box size, and especially near cut-off this 
increases the number of iterations required to solve the 
eigenvalue problem to a given accuracy (we used the well- 
known ARPACK solver to find the 11 lowest eigenmodes 
of the sparse matrix). We found that a box size of 6A 
to 7A gave the best compromise between accuracy and 
speed. Our non-uniform elements allowed us to have a 
high element density across the waveguide and in the 
trapping region, but a low density over the much larger 
box area, keeping the total number of degrees of freedom 
manageable. 

The fractional error e in the propagation constant k z 
was less than 1%, and the accuracy of the electric field 
strengths in the trapping region rs 3%, when we used 
N ~ 2000 elements. This was sufficiently accurate for the 
present work. Finding the bound modes of each waveg- 
uide parameter choice typically took between 3 and 20 
minutes of computing time on a Silicon Graphics R8000 
processor, depending on how close to cut-off the guide 
was. We tested the accuracy of the method by solving 
a cylindrical guide in an identical fashion with the same 
N and a very similar non-uniform grid, for which there 
are known field solutions |^8| . Figure || shows the prop- 
agation constant agrees with the analytics to within 1%, 
even close to cut-off. The convergence with N was mea- 
sured for the rectangular guide case, and found to be 
e ~ Af-T with 0.55 < 7 < 0.7. This is less than optimal 
for first-order elements (which have a maximum possi- 
ble convergence of 7 = 1), and is believed to be due to 
an inability of the bilinear functions to represent physi- 
cal in-plane E and H components at dielectric steps, or 
the weak field singularities which can physically occur at 
any exterior dielectric sharp edges (regardless of whether 
acute or obtuse) 

Future improvements to the method, which would in- 
crease the accuracy or the convergence rate 7, include 
using higher order elements (if done carefully, this could 
correctly represent physical E and H components at di- 
electric steps), and explicit modelling of the field singu- 
larities at guide corners using specialized elements. 



IV. FURTHER DECOHERENCE AND LOSS 
MECHANISMS 

A. Effects of surface interactions 

The EW trap has the benefit of creating high field 
gradients near a surface, but along with this comes the 
disadvantage that interactions with that surface that can 
alter the trap potential and even cause heating and loss 
of trapped atoms. 

An atom's change in potential near a surface is known 
as the van der Waals interaction (I <§; A) or the Casimir 
interaction (I A), depending on the distance I from 
the surface compared to A, the dominant wavelength re- 
sponsible for the polarizability of the atom (in our case 
of Cs this is the D line resonance, the same as our trap- 
ping resonance). There is a smooth cross-over from van 
der Waals (U ~ l~ 3 , which can be viewed as the atom's 
electrostatic interaction with the image of its own fluctu- 
ating dipole) to Casimir (U ~ Z~ 4 , which can be viewed 
as a retarded van der Waals attraction or equally well 
as an atomic level shift due to a cavity QED effect) at 
I « A/10 |}7). In the case of a perfect mirror surface, 
the full form is known for any I, but for a dielectric sur- 
face, the expression becomes much more complicated to 
evaluate j$8|. 

Since our trapping distances are larger than this cross- 
over point, we will use the Casimir form, which is correct 
for asymptotically large I, and is always an overestimate 
of the true potential M] . The dependence of the coeffi- 
cient with dielectric constant is complicated (3^^l[], but 
we will use the simpler approximate form given by Spruch 
and Tikochinsky p3), to give 



U Caa (l) 



1 



3 hca(0) 
8^ 47re Z 4 e+ (30/23)e 1 /2 + 7/23 



(4.1) 



(in the MKSA system) . This approximate form is known 
to be within 6% of the exact expression for any dielectric 
constant e |41|. Substituting the recently calculated |4C| ] 
static polarizability of cesium, a(0) = 399.9 a.u., gives 
a Casimir interaction coefficient of 4.9nK/im 4 for n g — 
1.56. 

Figure |l|b shows the effect of this potential on a typ- 
ical trap of depth 150 /iK and distance yo = 270 nm. It 
is clear that the change is negligible further than 100 nm 
from the surface, and a WKB tunneling calculation along 
this straight-line path (at x — 0) shows that even if all 
atoms that reach the surface stick, the loss rate from the 
first few transverse modes is entirely negligible. How- 
ever, care should be taken with the multi-mode regime, 
or in the case of high-p traps, since the tunneling via the 
corners of the "bean" shape may dominate for p > 0.4 
(Figure |). 

The issue of energy transfer to trapped atoms due to 
a finite (and possibly room) temperature nearby surface 
is far less well understood, and may be a problem with 
many surface-based particle traps, as discussed by Henkel 
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and Wilkens J42[ . However, since we are trapping neutral 
particles and the conductivity of our surface is low, we 
expect a decoherence rate negligible compared to that al- 
ready present from spontaneous absorption and emission 
cycles. 



B. Bending the waveguide 

It would be very useful to be able to carry atom beams 
along curved paths, by bending our atom waveguide in 
the plane of the substrate, without significant atom loss. 
Here we briefly estimate three limitations on the waveg- 
uide minimum bending radius (in decreasing order of le- 
niency): 1) the limit imposed by optical radiation leak- 
age, 2) the limit needed for incoherent atom transport, 
and 3) the limit needed for coherent atom transport in 
the transverse groundstate. This will give us an idea of 
the practicality of curved atomic guides. 

Firstly, whenever a dielectric optical guide has curva- 
ture, there is a loss rate (exponentially small in the cur- 
vature radius R), which can be viewed as tunneling out 
of the guide's 'potential well' induced by the addition of 
an effective centrifugal potential. In the limit R 3> W, 
the effective potential is linear with x (the radial coordi- 
nate), and the fractional loss per radian of curvature can 
be estimated , for instance using the one-dimensional 
WKB formula, to be 



C — exp 
A 



1 X 2 R 
6tt 2 LI 



(4.2) 



where L x is the evanescent decay length in the radial 
direction, and C is a constant of order unity. There- 
fore for negligible light loss at a 7r/2 bend we need 
R > 60ir 2 L x /\ 2 , typically a couple of tens of microns. 
This is so small chiefly because we are using an optical 
guide with a large refractive index step p7[ . 

Secondly, we consider atom loss from an incoherent 
beam with a transverse temperature ksT±_ (assumed 
small compared to the trap depth magnitude J7 max ), and 
a longitudinal kinetic energy E\\ . We call the approxi- 
mate spatial extent of the trap potential in the x direction 
2£, and restrict ourselves to one-dimensional classical mo- 
tion in this direction. When in a region of radius of curva- 
ture R ^> £, an effective centrifugal term adds to the trap- 
ping potential giving U{x) — Udi P (x) — 2E\\{x/R). This 
causes the atoms to 'slosh' towards positive x, only ever 
returning if there exists a point where U{x) > — i7 max 
for x > 0. We can estimate that this will happen 
if ^max > 2En(£/R), giving our lower limit on R as 
2£2£||/?7 max . In our design £ « 0.5 /mi, so if we choose 
R = 1 mm we can expect loss-free transport of a beam 
at a longitudinal kinetic energy up to 10 3 times the trap 
depth. 

Thirdly, to model coherent matter-wave propagation 
along a curved guide, we consider the amplitude for re- 
maining in the transverse ground-state, having passed 



into a curved section and back into a straight section. If 
again we assume one-dimensional x motion, and assume 
a harmonic potential U(x) — ^Mcu 2 (x — xq) 2 around the 
trap minimum, then the effect of curvature is to shift the 
minimum position from xq = to xq = 2E^/Mto 2 R. If 
this shift is much less than the characteristic ground-state 
size {Ti/Mujx) 1 ! 2 then the projection at each transition 
will be high, resulting in high flux transmission coeffi- 
cient. This gives R > 2E\\/ (hMw^) 1 / 2 as our condition, 
which for Eu/ks = lOmK (that is, Uii = 1.1ms ) and 
(v x /2ir = 40 kHz corresponds to R 0.45 mm. This limit 
is very conservative since we have not yet made use of the 
adiabatic condition fl <C uj x (where Q = v\\/R is the rate 
of change of direction of the guided atom), to design a 
waveguide path without discontinuities in the curvature. 

In conclusion, we have shown that it is possible to 
bend atoms both incoherently and coherently through 
large angles on a compact substrate structure of a few 
millimeters in size. 



V. CONCLUSION 

We have presented a novel substrate-based neutral 
atom waveguide, based on the optical dipole force, which 
combines the features of a planar far-detuned two-color 
evanescent trap jl3] with the ability to confine strongly 
along two axes. We utilized the differing vertical de- 
cay lengths of the two bound-mode polarizations of a 
submicron-sized optical waveguide near cut-off. We have 
shown that only a few milliwatts of guided laser power 
can create atomic potential depths ~ 100 fiK with trans- 
verse oscillation frequencies ~ 100 kHz, a coherence time 
~ 10 ms and a trap minimum 200-400 nm above the opti- 
cal guide surface, for Cs atoms in the mp = state. Laser 
powers greater than ten milliwatts can give a transverse 
mode spacing greater than the temperature of a Cs MOT, 
opening up the single-mode waveguide regime. The ad- 
vantages of guiding optical trapping fields on a substrate 
include mechanical stability and reliability, mass produc- 
tion and the potential for transport along complicated 
paths. 

We have given some design criteria for guided- 
lightwave two-color atom waveguides (chiefly the max- 
imization of the evanescent decay lengths, and of their 
normalized difference ol£), and shown that a substrate of 
low refractive index can be very beneficial. We modelled 
in detail the trapping potentials for a general rectangu- 
lar guide of index 1.56 above a unity- index substrate, and 
have shown that a realistic substrate choice of index 1.32 
poses few problems to the viability of the device. We pre- 
dict that the effect of the surface interaction is generally 
small, and that coherent guiding is possible around cor- 
ners of radii ~ 1 mm for a longitudinal velocity ~ 1 ms . 
We also believe that the magnetic part of the potential 
felt by nonzero atoms could be used to increase the 
depth and oscillation frequencies further. 







This preliminary work (specifically equation ( |2.3| )) in- 
dicates that utilizing detunings much larger than the 
15 nm we limit ouselves to here will be very advantageous 
for coherent guiding. We have only scratched the surface 
of the design variations possible; for instance, equaliz- 
ing the horizontal and vertical oscillation frequencies is 
yet to be attempted. The use of two polarizations is our 
solution to the problem of maximizing when the de- 
tuning is very small compared to the wavelength, but we 
suspect that there will exist other fruitful schemes where 
these are comparable (A ~ lu , very far-detuning) and 
where a large oll is caused simply by the different opti- 
cal cut-off conditions at the two wavelengths. We have 
reserved investigation of cooling schemes for future work 
(although this has already been demonstrated in an EW 
mirror jwj and proposed in EW traps p5| ). We believe 
that the potential shapes capable of being produced by 
guided waves on a substrate also include the possibility 
of funnel-type loading schemes and coherent atom cou- 
plers, allowing for a complete "integrated" atom-optical 
experiment on a substrate. 
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FIG. 1. (a) Shows trap geometry, dielectric guide dimen- 
sions, incoming laser polarizations, and the cartesian axes; (b) 
shows the trapping potential above the dielectric along a ver- 
tical slice at x = 0. The component due to red-detuned light 
(absolute value shown as dotted line) subtracts from that of 
blue-detuned light (dashed line) to give the total dipole po- 
tential C/dip (thin solid line). This is modified by the Casimir 
surface interaction (Section IV A), giving the final potential 
(thick line). Here the trap depth of 150 [iK and coherence time 
of 12 ms is generated in our design by 2 mW total guided laser 
power detuned by ±15 nm from the cesium D2 line. 
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FIG. 2. Numerically solved cut-off curves for a dielectric 
waveguide of n g — 1.56 (with n s = 1) as a function of its 
width and height (dashed curves) and, on the same axes, 
contours of the maximum trapping potential depth achiev- 
able at fixed total laser power (thin solid curves in rectangu- 
lar overlayed box region). Also shown is the symmetry line 
W — H (thin dash-dotted line). Note that the contours show 
depth increasing as H decreases, almost independent of W, 
near the suggested operating dimensions (shown as a solid 
ellipse). Cut-off is defined as reaching an effective refractive 
index n c g = k z /k — 1.05, except for i?fi and E\ x (thick 
dashed lines) which we show cut-off at n e g = 1.02. 
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FIG. 4. Comparison of our optical guide bound-mode nu- 
merical implementation against known analytic solutions, in 
the case of a free-standing dielectric cylinder of n g — 1.56. We 
used discretization and box-size identical to the rectangular 
guide case, and observe typical errors of ±0.5% in propagation 
constant for the first two modes. The mode naming conven- 
tion and analytic calculation follow Snitzer [glj; an asterisk 
indicates a doubly-degenerate mode. 



x (nm) 



1000 
800 
600 
400 



- 




dielectric 



x (nm) 



-1000 -500 500 1000 

FIG. 3. Trapping potential shapes in the xy plane, with 
guide dimensions of 0.83 /im by 0.21 /im and an index of 1.56. 
In each case the maximum depth is 100 /iK and the coherence 
time of atoms in the ground state is 19 ms, achieved with 
detuning ±15 nm from the cesium D2 line. The outer con- 
tour shows a depth of 3 pK, the cesium MOT temperature. 
Subsequent contours are spaced by 6 fiK. The plots illustrate 
the range of trapping distances acheivable: the upper trap 
(p = 0.4, using a total guided power of 1 mW) has a mini- 
mum 0.24 /jm from the surface; the lower trap (p = 0.2, total 
guided power of 7.5 mW) has a minimum 0.44 /im from the 
surface. 
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